% 计算各阶差商
clear; clc;
x = [-2,-1,1,2];
y = [5,3,17,21];
[p,q] = d_d(x,y);
